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Abstract 

Recently in the context of studies of the intrinsic charm content of the nu- 
cleon and of the rj' meson two groups have arrived at different results for the 
1/m 3 term of the heavy quark expansion for operator differing by the fac- 

tor of six. We show that the form of both results violates certain general con- 
ditions. Using the expression for the axial anomaly with the finite Pauli-Villars 
regularization we obtain a new expression for 1/m 3 term of the heavy quark ex- 
pansion for Q^Q. With this new result we obtain an estimate for the constant 
4' C) « -vh^Gl^W) « -2 MeV 

1. The analysis of the charm content in noncharmed hadrons often appeals to the heavy 
quark expansion of operators like Q75C} |], [|, [§. One starts with the calculation of the 
expectation value of (Q^Q) in the background gluon field A with the subsequent heavy 
quark expansion in inverse powers of m (we present the euclidean version) 

Recently two groups and JJ] came up with different results for the coefficient c 2 differing 
by the factor of six. Specifically the result of fl| for the large m expansion ([]]) is 

whereas the authors of obtain 

However, a simple argument shows that both results are not consistent with the general 
properties of the axial anomaly. Indeed, the diagonal matrix element of the operator in 
the lhs of eqs. (0), ([3]) is nothing else but the anomalous divergence of the axial current 
of a massless fermion with the heavy m fermion considered as a Pauli-Villars regulator 
(see e.g. the review [[§]) 



(d^tix)) = 2mSp(x| 75 |x) . (4) 

Eq. (|4]) is exact in the Pauli-Villars regular izat ion and should be valid in any order 
of the large m expansion. On other hand, one can write the local large m expansion for 
the axial current in the background gluon field 

i»> = E^M- ( 5 ) 

This means that in the large m expansion 

Sp(x| - * 75 b) = Y —d.X^ix) . (6) 

However, the term fabcF^F^F^ in the rhs of the large m expansions (Q), (|3|) can- 
not be represented as a total derivative of a local expression. Indeed, a straightforward 
calculation e.g. for the instanton field shows that 

d'xf abc F a P b ua F c atl £ . (7) 



A simple dimensional analysis shows that this nonvanishing contribution cannot come 
from the surface terms if the instanton field is taken in the regular gauge. Hence the 
integrand cannot be a total derivative. 

We have recomputed the 1/m 3 term of the expansion ([!]) and obtained the following 
result 

2mSp(x| 1 75 |x) = ^Sp color (f Xu F Xu ) + \ , (8) 

7^V M (v4) + m 87r^ v ' 9bir z m z 

Rfi = ^SPcoior (FxjsFx^ - 4Sp color ([V Q , F ua \F^ . (9) 

Here the fundamental representation of the gauge group is implied F\ v = F^ u ^-. This 
expression differs from both earlier results (0), (|3|) and explicitly has the form of a total 
derivative in agreement with the general form of the 1/m expansion (|J). In the case of 
the instanton field satisfying the self-duality equation F Xv = F Xu , we find 

^instanton = <^ S Pcolor {F Xu F Xu ) = , (10) 

where p is the size of the instanton. 

In principle the result (|J) can be also extracted from the calculation of the axial 
anomaly in the so-called finite- mode regularization in |J. Note that in that paper the 
parameter m played the role of the infrared regularization. Taking the ultraviolet cutoff 
of the finite mode regularization to infinity one obtains from eq. (5.7) in || the following 
expression 

2mSp(x| 1 75 1 x) = -^e^ Ap Sp color (F^F Ap ) + \ t^ uXp 

7 ^V M (A) + m \b-n z 9bix z m z 
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x Sp color {2[V a , [V a , F^Fxp + [V Q , F^] [V a , F Ap ] - 4iF Xfi F„ a F ap } . (11) 



After some simple algebraic transformations using the Bianchi identity, expression (|TT| ) 
can be brought to the form (H). 

2. Using the result (|8]) we can easily estimate the matrix element: 

(0|c7 M 75c|r ? '(g))=^4 c) g M . (12) 
Combining the anomaly equation 



Vcry^c = 2im c c l5 c + ^G%G^ a , (13) 



Air' 

and our result (rotated to the Minkowski space) we obtain the folowing estimate of 

f(c) 



the constant in the order m~ 

$ = -^M^^'M ■ ( 14 ) 



In this equation we use the "perturbative" normalization for the gluon field strength 
G a = F^/g and the term proportional to [V a ,G ua ], which vanishes in pure Yang-Mills 
theory, is neglected. The omitted term with the help of QCD equation of motion can be 
related to the following matrix element 

(y _ _ 

7Z<0|<? E ^fl^fW). (15) 

f=u,d,s 

A method of computation of such matrix elements in the framework of the instanton model 
of the QCD vacuum was developed in ref. 0. This method has already been applied to 
calculations of the nucleon matrix elements of certain mixed quark-gluon operators related 
to higher-twist corrections to deep-inelastic scattering 0. Rough order of maginitude 
estimate using the results of ref. J!]] indicates that the omitted term contributes (by 
absolute value) at the level of 0.3 MeV to the constant f$ . Full quantitative estimates 



l G%G^\ V ') = 0.056 GeV 3 , (16) 



of the matrix element (15) will be given elsewhere. 
Using the value 

<■ 

obtained in f8|, we get from ([14]) the following estimate: 

w -2 MeV . (17) 

Analogous estimate for the constant gives ~ —0.7 MeV, if one uses the value 
{®\%G a G a \ri) = 0.020 GeV 3 obtained in ref. gj. Let us note however that in the case 
of 77 the contribution of the omitted term can be of the same order, hence this estimate 
should be considered as pure. The full quantitative estimate of and in the order 
l/m| of heavy quark mass expansion will be given elsewhere. 

To summarize, our result for the m~ 3 term of the large m expansion (|XJ) agrees with 
the general properties of the axial anomaly in contrast to the expressions derived in [|[ 
and [Q] . It should be noted that the expressions of both M and were used by a number 
of authors for the analysis of the charm content of noncharmed hadrons (see e.g. Jl|, 0, [|). 
The structure of our result (^|) essentially differs from that of Jl] and which means that 
the results of papers relying on M and [fj must be reconsidered. 
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